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Abstract. We compute statistical equilibrium states of rotating self-gravitating fermions by maximizing the Fcrmi- 
Dirac entropy at fixed mass, energy and angular momentum. We describe the phase transition from a gaseous phase 
to a condensed phase (corresponding to white dwarfs, neutron stars or fermion balls in dark matter models) as we 
vary energy and temperature. We increase the rotation up to the Keplerian limit and describe the flattening of the 
configuration until mass shedding occurs. At the maximum rotation, the system develops a cusp at the equator. 
We draw the equilibrium phase diagram of the rotating self-gravitating Fermi gas and discuss the structure of 
the caloric curve as a function of degeneracy parameter (or system size) and angular velocity. We argue that 
systems described by the Fermi-Dirac distribution in phase space do not bifurcate to non-axisymmetric structures 
when rotation is increased, in continuity with the case of polytropes with index n > 0.808 (the Fermi gas at 
T = corresponds to n = 3/2). This differs from the study of Votyakov et al. (2002) who consider a Fermi-Dirac 
distribution in configuration space appropriate to stellar formation and find "double star" structures (their model 
at T — corresponds to n = 0). We also consider the case of classical objects described by the Boltzmann 
entropy and discuss the influence of rotation on the onset of gravothermal catastrophe (for globular clusters) and 
isothermal collapse (for molecular clouds). On general grounds, we complete previous investigations concerning 
the nature of phase transitions in self-gravitating systems. We emphasize the inequivalence of statistical ensembles 
regarding the formation of binaries (or low-mass condensates) in the microcanonical ensemble (MCE) and Dirac 
peaks (or massive condensates) in the canonical ensemble (CE). We also describe an hysteretic cycle between the 
gaseous phase and the condensed phase that are connected by a "collapse" or an "explosion" . 
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1. Introduction 



It is striking to observe that self-gravitating systems fol- 
low a kind of organization despite the diversity of their ini- 
tial conditions and their environment. This organization 
is illustrated by morphological classification schemes such 
as the Hubble sequence for galaxies and by simple rules 
which govern the structure of individual self-gravitating 
systems. For example, elliptical galaxies display a quasi- 
universal luminosity profile described by de Vaucouleur's 
i? 1 / 4 law and most of globular clusters are well-fitted by 
the Michic-King model. On the other hand, the rotation 
curves of spiral galaxies appear to be flat and this striking 
observation can be explained by the presence of a dark 
matter halo with a density profile decreasing as r~ 2 at 
large distances (Binney & Tremaine 1987). 



The question that naturally emerges is what de- 
termines the particular configuration to which a self- 
gravitating system settles. It is possible that their ac- 
tual configuration crucially depends on the conditions that 
prevail at their birth and on the details of their evolu- 
tion. However, in view of their apparent regularity, it is 
tempting to investigate whether their organization can 
be favoured by some fundamental physical principles like 
those of thermodynamics and statistical physics. We ask 
therefore if the actual states of self-gravitating systems 
are not simply more probable than any other possible con- 
figuration, i.e. if they cannot be considered as maximum 
entropy states (Chavanis 2002a). 

This thermodynamical approach may be particularly 
relevant for globular clusters and elliptical galaxies whose 
distribution function is close to the isothermal one (at 
least in the inner region of the system). In the case of 
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globular clusters, the relaxation proceeds via two-body 
encounters and this "collisional" evolution is governed by 
the Fokker-Planck-Landau kinetic equation for which a 
H-theorem is available (see Kandrup 1981 for a precise 
discussion). This collisional evolution can lead to the es- 
tablishment of a statistical equilibrium state. This equilib- 
rium state is usually described by the Michie-King model, 
a truncated isothermal, which takes into account the es- 
cape of high energy stars due to tidal forces (Michie 1963, 
King 1966). In fact, such equilibria are only metastable 
(Antonov 1962) and the evolution continues for longer 
times with the formation of binaries (Hcnon 1961). By 
contrast, for elliptical galaxies, two-body encounters are 
completely negligible (the corresponding relaxation time 
t co u exceeds the age of the universe by many orders 
of magnitude) and the galaxy dynamics is described by 
the Vlasov equation, i.e. collisionless Boltzmann equation. 
Yet, collisionless stellar systems can reach a metaequilib- 
rium state as a result of a violent relaxation (Lynden-Bell 
1967). The statistical prediction of this metaequilibrium 
state is complicated and depends on the initial condi- 
tions (Lynden-Bell 1967, Chavanis 2003a). Furthermore, 
the metaequilibrium state actually reached by the sys- 
tem may differ from the statistical prediction due to in- 
complete relaxation. Therefore, iJ-functions must be con- 
structed to account for the structure of elliptical galaxies 
(Tremaine, Henon & Lynden-Bell 1987, Hjorth & Madsen 
1993, Chavanis 2003a,c). 

The thermodynamics of self-gravitating systems 
started with the discoveries by Emdcn (1907), Bonnor 
(1956) and Antonov (1962) that a self-gravitating isother- 
mal gas can be in hydrostatic equilibrium only in a limited 
range of thcrmodynamical parameters (energy, tempera- 
ture and external pressure). Outside this range, the ther- 
modynamical potential (entropy, free energy or Gibbs en- 
ergy) has no extremum and the system is expected to col- 
lapse under its own gravity (Lynden-Bell & Wood 1968). 
A "core-halo" phase halting the collapse can be evidenced 
if we introduce appropriate small-scale cut-offs. Different 
regularizations have been proposed such as quantum de- 
generacy, hard spheres, soften potential,... Then, the shape 
of the caloric curve T(E) depends on an additional param- 
eter n related to the inverse of the small-scale cut-off (see 
Chavanis 2002b). For p — > +oo, one recovers the classical 
results of Emden, Bonnor and Antonov. As the cut-off in- 
creases, it is possible to describe phase transitions between 
a "gaseous" phase (independent on the small-scale cut-off) 
and a "condensed" phase with a core-halo structure. It is 
important to stress that statistical ensembles are not in- 
terchangeable for self-gravitating systems (Thirring 1970, 
Lynden-Bell & Lynden-Bell 1977, Padmanabhan 1990) 
and that phase transitions occur not only in the canon- 
ical ensemble but also in the microcanonical ensemble if 
the cut-off is sufficiently small (Stahl ct al. 1994, Chavanis 
2002b, Chavanis & Ispolatov 2002). Microcanonical phase 
transitions are a rather new phenomenon associated with 
the existence of metastable states (local entropy maxima) . 



In the previous studies, the system is supposed to 
be non-rotating. The thermodynamics of rotating stellar 
systems was studied by Lagoute & Longaretti (1996) in 
the framework of Michie-King models describing globu- 
lar clusters. The case of a slowly rotating isothermal gas 
confined within a spherical box was treated by Chavanis 
(2002c) using analytical methods introduced by Milne 
(1923) and Chandrasekhar (1933) for distorted poly- 
tropes. Recently, Fliegans & Gross (2002) and Votyakov 
et al. (2002) have considered the statistical mechanics of 
a rapidly rotating self-gravitating gas in relation with the 
process of star formation. They showed that phase transi- 
tions from axisymmetric to non-axisymmetric structures 
("double stars") occur at sufficiently high values of an- 
gular momentum. These authors provide an interesting 
phase diagram of rotating self-gravitating systems for ar- 
bitrary values of angular momentum and energy. They 
treat the particles as hard spheres (of size a) and con- 
sider a relatively large value of the small-scale cut-off. 
Fliegans & Gross (2002) evaluate the density of state 
g(E,L) numerically for a finite number N of particles 
in two-dimensions. Votyakov et al. (2002) use a mean- 
field approximation N » 1 and adopt a distribution 
/?($) corresponding to a Fermi-Dirac statistics in configu- 
ration space. This model is expected to take into account 
hard sphere effects, where the hard spheres are introduced 
heuristically to fix an upper limit to the local density. 

In this paper, we shall study the structure of rotating 
self-gravitating systems described by a Fermi-Dirac distri- 
bution function / = /(e) in phase space. This is the usual 
distribution function of quantum particles (fermions) due 
to Pauli's exclusion principle. This distribution function 
describes electrons in white dwarf stars (Chandrasekhar 
1942), neutrons in neutron stars (Hertel & Thirring 1971) 
and massive neutrinos in dark matter models (Bilic & 
Viollier 1997). Incidentally, this distribution function also 
occurs in the context of the violent relaxation of collision- 
less stellar systems (Chavanis & Sommeria 1998). In the 
latter, the statistical approach refers to a coarse-graining 
of the distribution function and the effective exclusion 
principle is a consequence of the Liouville theorem associ- 
ated with the collisionless nature of the evolution (Lynden- 
Bell 1967, Chavanis 2002d). 

The distinction between the distribution used by 
Votyakov et al. (2002) and the ordinary Fermi-Dirac dis- 
tribution function in phase space is crucial, regarding the 
structure of the condensed phase, when the system is ro- 
tating. At T = (or E = E min ), the former leads to 
a homogeneous body (p = Cst.) while the latter leads 
to a polytrope with index n = 3/2 corresponding to a 
classical white dwarf star (p = Kp 5 / 3 ). Now, it is well- 
known that a uniformly rotating homogeneous mass is 
axisymmetric (spheroidal) for slow rotations (Maclaurin 
sequence) but bifurcates to non-axisymmetric configura- 
tions at sufficiently high rotations (Chandrasekhar 1969). 
These non-axisymmetric structures are usually ellipsoidal 
(Jacobi sequence) but more complex structures (pear- 
shaped, dumbbell,...) have also been found. The fission of 
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the self-gravitating fluid into a binary state is also possible 
(Hachisu & Eriguchi 1984) and these binary states appear 
in the study of Votyakov et al. (2002) at low temperatures 
and high rotations. On general grounds, the fission insta- 
bility has been extensively discussed in relation with the 
formation of binary stars and the formation of the moon 
(in the non-symmetrical case). However, the fission sce- 
nario may not be the relevant mechanism and it is almost 
abandoned at present. Binary stars are more likely to re- 
sult from the fragmentation of the molecular cloud (Boss 
1988) and the moon has probably been detached from 
the protoearth after a collision with a Mars-sized body 
(Boss 1986). In contrast to incompressible fluid masses, a 
polytropic configuration with index n > 0.808 does not 
bifurcate to non-axisymmetric structures when angular 
velocity is increased; it remains axisymmctric until the 
Keplerian limit is reached (James 1964). At that point 
there is mass shedding at the equator (with the formation 
of a cusp) and disruption of the system. This concerns 
in particular polytropes of index n — 3/2 (white dwarfs) 
which are the zero temperature limit of the Fermi distri- 
bution (Fowler 1926). This picture probably persists at 
finite temperature where we have, typically, a Fermi con- 
densate surrounded by a "vapour". Therefore, the phase 
diagram corresponding to the Fermi-Dirac statistics in 
phase space is expected to differ from the one calculated 
by Votyakov et al. (2002). In this paper, we propose to 
compute this phase diagram for arbitrary rotation and de- 
generacy parameter (or system size) in both microcanon- 
ical and canonical ensembles. We shall restrict ourselves 
to axisymmetric single-cluster configurations because our 
code is limited to such configurations. According to James' 
result, if other structures exist (double-clusters, rings,...), 
they should bifurcate dis continuously from the spheroidal 
sequence as rotation increases. Such bifurcations will not 
be considered in this paper. The present study com- 
pletes the description of non-rotating or slowly rotating 
fermionic structures discussed by Chavanis (2002b, c). 

2. Statistical mechanics of rotating 
self-gravitating systems 

2.1. The maximum entropy state 

Consider a system of N particles, each of mass m, 
interacting via Newtonian gravity. We allow the sys- 
tem to have a non-vanishing angular momentum. Let 
f(r,v,t) denote the distribution function of the system, 
i.e. f(r,v,t)d 3 rd 3 v gives the mass of particles whose po- 
sition and velocity are in the cell (r, v; r + d 3 r, v + d 3 v) at 
time t. The integral of / over the velocity determines the 
spatial density 



(i) 



The mass and angular momentum of the configuration are 
given by 



M 



-J 



pd 3 r, 



(2) 



L = J f rxv d 3 rd 3 v. 



(3) 



On the other hand, in the mean-field approximation, the 
energy can be expressed as 



E = ^ I fv 2 d 3 rd 3 v + i / p$d 3 r = K + \V. 



(4) 



where K is the kinetic energy and W the potential energy. 
The gravitational potential $ is related to the density by 
the Newton-Poisson equation 



A$ = AirGp. 



(5) 



Let us consider different systems of astrophysical in- 
terest that can be approached by statistical mechanics, (i) 
The equilibrium states of collisional stellar systems (e.g., 
globular clusters), reached for t > t re i ax <~ ^n^to (where 
to is the dynamical time), are obtained by maximizing the 
Boltzmann entropy (Antonov 1962, Ogorodnikov 1965, 
Lynden-Bell & Wood 1968) 



/ 



/ln/d 3 rd 3 v, 



(6) 



at fixed mass M, energy E and angular momentum L. 
This is a condition of thermodynamical stability for the 
Hamiltonian N-stars system. For globular systems, only 
the microcanonical ensemble makes sense and equilibrium 
states exist only above the Antonov energy. Furthermore, 
these equilibria are only metastable (but long-lived) since 
the Boltzmann entropy has no global maximum, (ii) 
Collisionless stellar systems (e.g. elliptical galaxies) can 
reach a metaequilibrium state on a very short time scale 
~ tu as a result of a violent relaxation. This metaequi- 
librium state is a stationary solution of the Vlasov equa- 
tion. Lynden-Bell (1967) has tried to predict this metae- 
quilibrium state by resorting to a new type of statistical 
mechanics. Unfortunately, his prediction is limited by the 
problem of incomplete relaxation (e.g., Hjorth & Madscn 
1993). Therefore, H-functions must be introduced to ac- 
count for the structure of galaxies (Tremaine, Hcnon & 
Lynden-Bell 1987, Chavanis 2003a, c). It can be shown that 
the maximization of a H-function 



rd 3 v, 



(7) 



at fixed M, E and L, where C(f) is a convex func- 
tion (i.e. C"(f) > 0) determines nonlinearly dynamically 
stable stationary solutions of the Vlasov equation (Ipser 
& Horwitz 1979). This is similar to a condition of mi- 
crocanonical stability in thermodynamics where S plays 
the role of a generalized entropy (Chavanis 2003a). Note 
that according to the Jeans theorem (Binney & Tremaine 
1987), other solutions of the Vlasov equation can be con- 
structed which do not satisfy this criterion, (iii) The equi- 
librium properties of stars are obtained by coupling the 
hydrostatic equations to a local equation of state. In sim- 
plest models, it is assumed that the star is barotropic, i.e. 
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p = p{p). We introduce the energy functional 



1 



o P 

2 



p^-d 3 r — f2 j p r x u d 3 r, 



(8) 



where the first term is the internal energy, the second the 
gravitational energy, the third the kinetic energy of the 
bulk motion and the third the energy of rotation. It can 
be shown that the minimization of W[p] at fixed mass M 
and angular velocity H determines nonlinearly dynami- 
cally stable stationary solutions of the Euler- Jeans equa- 
tions. This is similar to a condition of canonical stability in 
thermodynamics where VV plays the role of a generalized 
free energy (Chavanis 2003a). (iv) Finally, the interstellar 
medium can be considered as a gas in thermal equilib- 
rium with a radiation background (de Vega & Sanchez 
2002). Therefore, for this system, the canonical descrip- 
tion makes sense. The stability of the gas against grav- 
itational collapse can be studied thcrmodynamically by 
minimizing the free energy F = E — TS associated with 
the Boltzmann entropy © or dynamically by using the 
Jeans instability criterion (see Chavanis 2002e). This is 
the traditional approach to understand star formation. 

In this paper, we shall consider a system of self- 
gravitating fermions. They can be electrons in white dwarf 
stars, neutrons in neutron stars or massive neutrinos in 
dark matter models. Their statistical equilibrium state is 
obtained by maximizing the Fermi-Dirac entropy 



S = 



J_ 



In 1 



'/o 



d 6 rd 6 v, (9) 



at fixed mass, energy and angular momentum. In this for- 
mula t]q = (2s + l)m 4 / (2irft,) 3 is the maximum allowable 
value of the distribution function fixed by Pauli's exclu- 
sion principle (s is the spin). For the quantum gas, we 
shall consider either microcanonical or canonical equilib- 
ria. Incidentally, the Fermi-Dirac entropy © also occurs 
in the statistical theory of violent relaxation for collision- 
less stellar systems developed by Lynden-Bell (1967). It 
corresponds to the two-levels approximation /o € {0, 770} 
of his theory, where /o is the initial distribution func- 
tion (Chavanis & Sommeria 1998). In the non degenerate 
limit / <C rjo, the Fermi-Dirac distribution (|5J) returns the 
Boltzmann entropy © as a special case. Therefore, our 
study can apply either to quantum particles, stellar sys- 
tems or to the interstellar medium. Of course, each system 
corresponds to a different limit of the theory characterized 
by a different degree of degeneracy and rotation. 

Looking for critical points of entropy 10 at fixed mass, 
energy and angular momentum and introducing Lagrange 
multipliers a, (3 and — /3fi for each constraint, we obtain 
the Fermi-Dirac distribution function 



/ 



Vo 



1 + Ae^-^j) ' 



(10) 



where e = \ + <I> and j = r x v are the energy and the an- 
gular momentum of a particle by unit of mass and (3 = 1/T 



is the inverse temperature (we include m or rjo in the def- 
inition of T). The same distribution IjlOfl is obtained by 
extrcmizing the free energy J = S — /3S+/3J7-L at fixed (3, 
f2 and M. However, the second order variations of S and J 
(under respective constraints) can have a different sign im- 
plying an inequivalence of statistical ensembles (Chavanis 

2003c) . Introducing the Jacobi energy e j = e — f2 • j and 

2 

noting that e j — ^ + / , where w = v — $7 x r is the 
relative velocity and & e ff = 3> — ij(fixr) 2 is the effective 
potential accounting for inertial forces, we can rewrite Eq. 
(|10fl in the form 



/ = 



'/(> 



1 + Ae^T" 



(11) 



We note that the "most probable" form of rotation is a 
rigid rotation f2. The gravitational potential $ is defined 
within an arbitrary gauge constant. We shall take the or- 
dinary convention $ — > at infinity. 

By integrating Eq. Ijll|l over the velocity, we find that 
the density is given by 



P : 



^3/2 J l/2l Ae )> 



where /„ is the Fermi integral 



In(t) = 



1 + te x 



-dx. 



(12) 



(13) 



The equilibrium configuration is then obtained by solving 
the Fermi-Poisson equation 



(14) 



and relating the Lagrange multipliers A, (3 and $7 to the 
constraints M, E and L. We stress that, for long-range 
systems, the mean- field approximation is exact in a suit- 
able thermodynamic limit (see Sec. 12.4(1 so that our "ther- 
modynamical" approach based on the maximization of 
thermodynamical potentials is rigorous and simpler than 
the "statistical" approach based on functional integration 
(Laliena 1999, Votyakovet al. 2002). The relation between 
the two points of view is discussed in Chavanis (2003a). 

2.2. The box model 

It is easy to see that Eq. (I14fl has no physical solution. 
Indeed, the density decreases as r~ 2 at large distances 
so that the total mass of the configuration is infinite in 
contradiction with our starting hypothesis. This means 
that there is no maximum entropy state in an unbounded 
domain. A self-gravitating system at non-zero temper- 
ature has the tendency to "evaporate" . Therefore, the 
statistical mechanics of self-gravitating systems (classical 
or quantum) is essentially an out-of-cquilibrium problem. 
However, evaporation is a slow process and, in practice, 
the system is found in a quasi- equilibrium state with a 
"core-halo" structure. The core of the system is approx- 
imately isothermal while the halo departs from isother- 
mality. Therefore, the distribution function (|10f> must be 
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modified at high energies. Some modified isothermal dis- 
tributions have been introduced in the case of globular 
clusters (Michie 1963, King 1966) subject to the tides of 
a nearby object (typically a spiral galaxy) and in the case 
of elliptical galaxies undergoing incomplete violent relax- 
ation (Hjorth & Madsen 1993). 

Another astrophysically plausible distribution function 
which takes into account Pauli's exclusion principle and 
tidal effects is given by (Chavanis 1998) 



/ = m- 



A + e-P* 



(15) 



for e < e m and / = for e > e m (where e m is the escape 
energy). This distribution function can be derived from 
a Fokker-Planck equation appropriate to self-gravitating 
fermions. It can describe dark matter halos limited in ex- 
tension by tidal effects (and, possibly, certain collisionless 
stellar systems). In the classical limit, it reduces to the 
Michie-King model. The distribution function (|15|l will be 
studied specifically in another paper. For the present, we 
shall keep the full Fermi-Dirac distribution function (jll|l 
and confine the system within a spherical box of radius 
R so as to avoid the infinite mass problem while preserv- 
ing the rotational symmetry of the system. Typically, the 
size of the box delimitates the physical extent of the sys- 
tem under consideration. Putting the system in a box is a 
simple and crude way of handling the infinite mass prob- 
lem. This is sufficient, in a first step, if we are just in- 
terested in describing phase transitions between gaseous 
and "core-halo" states. If, on the other hand, we wish to 
make more realistic models of self-gravitating systems, the 
spatial confinement must be given further consideration. 

2.3. Asymptotic limits T = and T -> +oo 

It will be convenient in the following to work with dimen- 
sionlcss variables defined by 

vp $ 



R' 



A 



ER 

GAP ' 

L' = 



(M/R 3 ) ' 
(3GM 



R 



L 



VGM 3 i?' 



/' 



(GM/R) ' 

\gm) ' 

7? \/5127r 4 G 3 Mi? 3 . (16) 



ST 



This is equivalent to setting M = R = G = 1 in the 
dimensional equations and we shall adopt this convention 
in the following. Then, the degeneracy parameter is given 
by (i — 16\/27r 2 ?7o- We also define tp = /3$ and u> = y/J3£l. 
Denoting by tpo the value of ip at r = and introducing 
the uniformizing variable k = Ae^° , the Fermi-Poisson 
equation (jl4|l becomes 



Aip = ,iVr/ 1/2 (fc^*^ u2s2 ), 



(17) 



where s — r sin 9 denotes the distance to the axis of ro- 
tation. Equation l|17l) must be solved with the condition 
tp — > at infinity. Before solving this equation in the gen- 
eral case, we first discuss important asymptotic limits. 



Forfc 
integral 

1. 



+oo, we can use the limiting form of the Fermi 



7 n (t) ~ -r(n 



1), (t _►+(»), 
and we obtain the Boltzmann-Poisson equation 



Atp 



2k 



(18) 



(19) 



for a rotating isothermal gas. This corresponds to the high 
temperature limit of the Fermi distribution. Indeed, for 
T ^> 1 the distribution function l|llf> is Maxwellian 



Alternatively, for k 
the Fermi integral 



(20) 

0, we can use the limiting form of 



I n {i) 



{-hit) 



n+l 



(t-0), 



n + l ' 

and obtain the Lane-Emden equation 
AV;= ^VT(-lnfc-^ + ^ + iw 2 S 2 ) 3/2 , 



(21) 



(22) 



for a rotating polytrope of index n — 3/2. This corre- 
sponds to the low temperature limit of the Fermi distri- 
bution in which the structure is completely degenerate 
(white dwarf). Indeed, for T = 0, the distribution func- 
tion 1)11(1 is a step function 



/ = VoH{ej - e F ), 



(23) 



where ej is the Jacobi energy, tp = — (l//3)lnA is the 
Fermi energy and H is the Heaviside function. Therefore, 
the Fermi-Dirac distribution function connects continu- 
ously isothermal and polytropic distributions for high and 
low temperatures. 

We recall that for the self-gravitating Fermi gas at T = 
(white dwarf), the density vanishes at a finite radius i?*. 
The mass-radius relation can be written (Chandrasekhar 
1942) 



MRi = 

with x 
to 



X 



(24) 



5.9723 1CT 3 . The condition R* < R is equivalent 

(25) 



fj, > Ai* = v / 512tt 4 x ^ 17.259 
More generally, we have 



R 



3/2 



(26) 



Therefore, the degeneracy parameter fi is related to the 
ratio of the system size R to the radius R„ of a white 
dwarf star with mass M. Accordingly, a small value of /x 
corresponds to a large "effective" cut-off (played by Pauli's 
exclusion principle) or, equivalently, to a small system size. 
Alternatively, a large value of /i corresponds to a small 
"effective" cut-off or a large system size. 
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2.4. The thermodynamic limit 

Before going further, we shall discuss the thermodynamic 
limit of self-gravitating systems. In order to precisely de- 
fine the thermodynamic limit, we must introduce a small- 
scale regularization otherwise the density of states g(E) 
diverges (Padmanabhan 1990). The most physical regular- 
ization is to consider quantum mechanics effects. Indeed, 
compact objects such as white dwarf stars and neutron 
stars owe their stability to quantum mechanics. The ther- 
modynamic limit corresponds to N ~ * +oo such that A = 
-ER/G N 2 m 2 , 7? = f3GNm/R, Q,' = Vl{R 3 /GNm) 1 ' 2 , 
V = L/VGN 3 m 3 R and n = (gm A /h 3 )^ 512tt 4 G 3 NmR 3 
are kept fixed. It is natural to fix m, h and G. Then, the 
strict thermodynamic limit is N — > +oo with 

R~N~ 1/3 , E~N 7/3 , /3~N- 4/3 , 
fl~N, L~N 4/3 . (27) 

This scaling was indicated in a previous paper (Chavanis 
2002b). When this limit is taken, S ~ N, J ~ N and the 
mean-field approximation is exact (note that the usual free 
energy F = E — TS = —J/fl scales as TV 7 / 3 ). The scaling 
MR 3 ~ 1 corresponds to the "white dwarf" mass-radius 
relation (|24p. The fact that R — > as N — > +oo accounts 
for the natural tendency of self-gravitating systems to col- 
lapse (note that when M reaches Chandrasekhar's limiting 
mass, Newtonian mechanics loses its relevance). We shall 
call (|27|l the quantum, or white dwarf, thermodynamical 
limit (QTL). It is particularly relevant for compact objects 
forming the condensed phase. 

If we now consider the classical limit h = 0, we see that 
/i becomes infinite so that it disappears from the problem. 
Then, the classical thermodynamic limit (CTL) is N — ► 
+00 with fixed A = -ER/GN 2 m 2 , r] = f3GNm/R, W = 
VL(R 3 /GNm) 1 / 2 , V = L/^/GN 3 m 3 R. We note that there 
is a freedom since we have to fix a parameter in order to 
determine the others. If we assume (3 ~ 1, the classical 
thermodynamic limit is N —* +oo with 

R ~ N, E ~ N, (3 ~ 1, 
n~N-\ L~N 2 . (28) 

This has been called the "dilute limit" by de Vega & 
Sanchez (2002). Again, S ~ N, J ~ N (and now F ~ N) 
so that the mean-field approximation is exact in that 
limit. We recall that for h = 0, the density of state di- 
verges due to the existence of collapsed configurations. 
However, for sufficiently large T and E/N, we know that 
there exists long-lived metastable gaseous states (local en- 
tropy maxima). For N ^> 1, these metastable states have 
a very long lifetime, increasing as exp(N), so they con- 
stitute the most relevant astrophysical structures (Katz 
& Okamoto 2000, Chavanis & Ispolatov 2002, Chavanis 
2003a). Globular clusters are in such metastable states. 
These gaseous states are independent on the small-scale 
cut-off so they are described adequately by the dilute limit 
Q28[). The consideration of these metastable states may 
solve the problems raised by Laliena (2003). Collapsed 



configurations (binaries in MCE and Dirac peaks in CE), 
which correspond to global maxima of entropy or free en- 
ergy (i.e., true equilibrium states) must be discarded "by 
hands" because they are reached for unphysically large 
times t — > +oo. They can be reached for accessible times 
only if the temperature and the energy are sufficiently 
small (T < T c or E < E c ) because, in that case, the 
gaseous phase ceases to exist and the system undergoes 
gravitational collapse. Clearly, the end-state of the col- 
lapse depends on the nature of the small-scale cut-off (e.g., 
quantum degeneracy) and we must consider the thermo- 
dynamical limit (|27|l if we want to describe these con- 
densed states. 

If instead of fermions we consider a hard sphere gas, 
the filling factor is [i = R/aN 1 / 3 where a is the hard 
sphere radius (see Chavanis 2002). Then, fixing a, m and 
G, the strict thermodynamic limit is N — > +oo with 

R-N 1 / 3 , E ~ N 5 / 3 , (3~N- 2 / 3 , 

f2 ~ 1, L ~ N 5/3 . (29) 

Finally, if we consider a soften potential with soften 
parameter r$, the cut-off parameter is /i = R/r$ (see 
Chavanis & Ispolatov 2002). Then, fixing ro, m and G, 
the strict thermodynamic limit is N — ► +oo with 

E~N 2 , f3~N-\ 
n~N 1/2 L~N 3/2 . (30) 

The dilute limit (|28|l is of course unchanged and it de- 
scribes gaseous metastable states independent on the 
small-scale cut-off. Note that other combinations keeping 
A, r), £1' and V fixed can be considered (e.g., G ~ 1/N for 
fixed m, R, /?, E/N, Cl and L/N orm~l/JV for fixed G, 
R, f3, E, £1 and L). Of course, the limits considered above 
are only formal. In reality, one does not have to make N, 
R,... go to +oo or 0. The precise values taken by these 
parameters are fixed by the physics of the system under 
consideration (molecular clouds, globular clusters, white 
dwarfs, neutron stars, dark matter halos,...). Note that 
even for small values of N, the mean-field approximation 
is accurate so that the thermodynamic limit makes sense 
even if we are not in the strict, mathematical, conditions 
of its application. 

2.5. Integral constraints 

According to the classical work of James (1964), rotat- 
ing polytropes with index n > 0.808 do not bifurcate to 
non-axisymmetric structures when rotation is increased. 
Instead, they develop a cusp at the equator. Since the 
Fermi gas at T — is equivalent to a polytrope of in- 
dex n = 3/2, it is not expected to bifurcate to non- 
axisymmetric structures (in contrast to homogeneous bod- 
ies). We expect this property to persist at non-zero tem- 
perature since, at T ^ 0, we have essentially a poly- 
trope 3/2 surrounded by a dilute halo. We do not ex- 
clude the possibility of discontinuous bifurcations to non- 
axisymmetric structures or rings. However, in this paper, 
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we shall restrict ourselves to single-cluster axisymmetric 
solutions of Eq. ifTTj) . 

In that case, the angular momentum can be written 
L = 7f2 where / is the moment of inertia 



I = 2n 



pr A {l — /i 2 )drdfi, 



(31) 



where fi = cos 9. The energy can be deduced from the 
Virial theorem (see, e.g., Chavanis 2002c) 



E= -K+ i> pr-dS. 



(32) 



The kinetic energy K incorporates a rotational contribu- 
tion 



K rot = 1l • n = ~/n 2 , 

and a thermal contribution 
Kth = \ f pd 3 r, 



(33) 



(34) 



where p = -g J fw 2 d 3 w is the local pressure. Using Eq. 
<|11|) . it can be written 



or, in dimcnsionless form, 

07T ' 



(35) 



(36) 



Equations l|12|) and (|35() determine the equation of state 
of the Fermi gas at non-zero temperature (in a parametric 
form). Since p = p(p), this gas is barotropic. 

2.6. A short description of the numerical procedure 

Using the variables introduced previously, the Gauss the- 
orem can be written as 



W^-dS 



47T 



(37) 



where the surface integral is computed on the spherical 
box. Now, using the decomposition of the gravitational 
potential in normalized spherical harmonics, i.e., 



l,m 



it turns out that 
4tt 



T 



^o(l) 



(38) 



(39) 



Therefore, Eq. (|17ll depends only on three parameters /x, 
and k. The degeneracy parameter \x and the angular 
velocity Q specify a series of equilibria parameterized by 
the uniformizing variable k. Given these three parame- 
ters, the Fermi-Poisson equation (|17ll can be solved by an 
iterative procedure. For that purpose, the radial compo- 
nents of the potential, i.e. ipi m (r), are discretized on the 



Gauss-Lobatto grid which ensures spectral convergence. 
Typically, the resolution used for the following numerical 
results was 100 radial grid points and a decomposition on 
16 spherical harmonics. 

The temperature T(/x, f2, k), the energy E(n, fi, k) and 
the angular momentum £(/i, fl, k) are determined by Eqs. 
(|39")l . and (j3T|) . It is then possible to compute the 
caloric curves E(T) for given [i and f2 (or L) by eliminating 
k from the previous relations. 



3. Caloric curves for low values of the degeneracy 
parameter (high cut-off/small system) 

We first describe the figures of equilibrium of the rotat- 
ing self-gravitating Fermi gas for a degeneracy parameter 
/i = 100. This value is smaller than the microcanonical 
critical point Hmtp — 2670 above which a first order mi- 
crocanonical phase transition appears (Chavanis 2002b). 
Hence, for /i = 100 (small system), the influence of the 
small-scale cut-off (played here by the exclusion princi- 
ple) is strong in average and stabilizes the system with 
respect to the gravothermal catastrophe (Lynden-Bell & 
Wood 1968). 

3.1. Non-rotating configurations 

The caloric curve of non-rotating self-gravitating fermions 
has been discussed by Chavanis (2002b) in both canonical 
and microcanonical ensembles and for arbitrary values of 
the degeneracy parameter /j, (or system size R). This com- 
pletes the work of Hertel & Thirring (1971) who worked in 
canonical ensemble and considered small system sizes. For 
non-rotating systems, the results obtained with fermions 
are similar to those obtained with other small-scale regu- 
larizations, the small-scale cut-off playing the same role as 
the inverse of the degeneracy parameter (Stahl et al. 1994, 
Follana & Laliena 2000, Chavanis & Ispolatov 2002). 

We first discuss the microcanonical ensemble corre- 
sponding to an isolated system with fixed energy. For 
fi < j-iMTP = 2670, there exists only one extremum of 
entropy for each energy and it is a global entropy max- 
imum (GEM). The caloric curve is plotted in Fig. ^ ( m 
this figure fi = 1000). For high energies, the density is 
almost homogeneous (Fig. [5^). This forms the "gaseous" 
phase. For energies close to the minimum accessible en- 
ergy E min = —6.42 10~ 2 // 2 / 3 (ground state), the temper- 
ature goes to zero and the system is completely degen- 
erate (Fig. |5^). It has the same structure as a classical 
white dwarf star, i.e. a polytrope of index n = 3/2. This 
forms the "condensed" phase. For intermediate energies 
(Econd < E < E gas ), corresponding to the region of neg- 
ative specific heats C = dE/dT < 0, the system has a 
"core-halo" structure with a degenerate nucleus (fermion 
ball) and a dilute isothermal atmosphere (Fig. 01) ■ This is 
like a solid condensate embedded in a vapor in usual phase 
transitions. This "mixed phase" is intermediate between 
the pure condensed phase and the pure gaseous phase. 
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Fig. 1. The caloric curve for a non-rotating system with 
/i = fOOO < Hmtp = 2670. fn the microcanonical ensem- 
ble (fixed E), the system slowly evolves from the gaseous 
phase to the condensed phase, passing by a mixed phase, 
as energy decreases. In the canonical ensemble (fixed T), 
the system undergoes a sharp phase transition at a crit- 
ical temperature (see text). Below T c , the gaseous states 
undergo a gravitational collapse (Jeans instability); above 
T*, the condensed states undergo an explosion. Between, 
T t and T c (resp. T„) the gaseous (resp. condensed) states 
are metastable. By varying the temperature of the system 
between T» and T c (spinodal points), we can generate an 
hysteretic cycle between the gaseous and the condensed 
phase. 



We now describe the canonical ensemble in which the 
temperature T is fixed instead of the energy E. Thus, we 
need to find maxima of free energy J = S — f3E at fixed 
(3 and M . We recall that, for an isolated self-gravitating 
system, the validity of using the canonical ensemble to de- 
scribe the statistics of a sub-system cannot be established 
from the microcanonical ensemble since energy is not ad- 
ditive (Padmanabhan 1990). Therefore, the notion of ther- 
mostat has to be defined carefully. The canonical ensem- 
ble can make sense for self-gravitating systems subject to 
additional short-range interactions. They can constitute 
a thermostat imposing a temperature. In that case, the 
usual canonical ensemble interpretation holds. For exam- 
ple, a self-gravitating system in contact with a radiation 
background (e.g., a molecular cloud) is usually treated 
in the canonical ensemble. This is the traditional way of 
studying stellar formation in relation with Jeans insta- 
bility. Note that there exists a precise theoretical model 
where the canonical ensemble is rigorously justified. This 
is the self-gravitating Brownian gas model (Chavanis et 
al. 2002, Sire & Chavanis 2002, 2003). In this model, the 
temperature has a precise physical interpretation since it 
is related to the strength of the stochastic force ( "kicks" ) 
acting on the particles. 



In the canonical ensemble, the region of negative spe- 
cific heat is unstable as it corresponds to saddle points 
(SP) of free energy (Katz 1978, Chavanis & Sommeria 
1998). It is replaced by a canonical first order phase tran- 
sition connecting the gaseous phase for T > T t (/x) to the 
condensed phase for T < T f (/i). Therefore, the mixed 
"core-halo" states with C < are thermodynamically for- 
bidden in the canonical ensemble. The gaseous states with 
temperature T c < T < Tt are metastable (local maxima 
of free energy LFEM) but they are long-lived. They can 
play an important role reminiscent of a supercooled state 
of the Van der Waals gas. At T < T c they undergo a grav- 
itational collapse (zeroth order phase transition) coincid- 
ing with Jeans instability criterion (Chavanis 2002e). For 
quantum particles, the outcome of this collapse is the for- 
mation of a "fermion ball" containing almost all the mass. 
In the classical limit [i — > +oo, the isothermal collapse 
leads to a Dirac peak (see Sec. EJ. The condensed states 
with temperature Tt < T < T* are also metastable and 
long-lived (if they are initially prepared in such states). 
They will undergo an "explosion" , reversed to the col- 
lapse, if they are heated above T*. This explosion trans- 
forms the dense core into a relatively uniform mass dis- 
tribution. By varying the temperature between T* and T c 
we can generate an hysteretic cycle in the canonical en- 
semble. This hysteretic cycle was described qualitatively 
in Chavanis (2003b). Recently, it has been followed nu- 
merically by using a model of self-gravitating Brownian 
fermions (Chavanis et al. 2003). The turning points of 
temperature where the metastable states terminate are 
called (canonical) spinodal points in the langage of phase 
transition. The first order canonical phase transition disa- 
pears for \x < ficTP — 83 (canonical critical point), i.e. for 
very high cut-off values. This is simply because the small- 
scale "repulsion" prevails over gravity. More details on 
phase transitions in self-gravitating systems can be found 
in Chavanis (2002b). 



3.2. Slowly rotating configurations 

The slowly rotating configurations of an isothermal gas 
(possibly degenerate) can be determined semi-analytically 
by using perturbative expansions of the Boltzmann- 
Poisson or Fermi-Poisson equations in terms of the di- 
mensionless rotation parameter v — f2 2 /2nGpo (Chavanis 
2002c). To lowest order in the expansion, the effect of ro- 
tation is just to flatten the system and to shift the onset 
of instability (see Sec. 03). For T = (or E = E min ), we re- 
cover the distorted polytrope of index n — 3/2 studied by 
Milne (1923) and Chandrasekhar (1933). In the region of 
negative specific heats, this "rotating fermion ball" is sur- 
rounded by a gaseous halo. For high energies, the system is 
non degenerate. Hence, the description of the equilibrium 
phase diagram is very similar to that of the non-rotating 
Fermi gas. 
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Fig. 2. Plot of the density (logarithmic scale) in a meridional section of the box. The solid line marks the level where 
density is 0.05, the central density being unity. The first row shows configurations at very low temperature/energy 
(condensed phase). The second row shows configurations in the region of negative specific heats (mixed phase). Such 
structures are forbidden (unstable) in the canonical ensemble. The third row shows configurations at high tempera- 
ture/energy (gaseous phase). Alternatively, the first column corresponds to non-rotating systems, the second column 
to slowly rotating systems and the third column to systems rotating close to the Keplerian limit. 



3.3. Spheroidal structures 

For higher values of rotation, the flattening increases and 
the system takes a spheroidal structure. For E — ► E m i n , 
we have a rapidly rotating "fermion spheroid" (Fig. I2>) . 
This structure is completely degenerate and highly con- 
centrated. In particular, it is insensitive to the confin- 
ing box. As energy increases further, the spheroid is only 
partially degenerate and extends to larger and larger dis- 
tances (Fig. Et) ■ The rotation affects in priority the low 
density external region while the dense core of the con- 
figuration remains almost spherical. In the region of neg- 
ative specific heats (E con d < E < E gas ), the "spheroid" 
is surrounded by a diffuse halo. This corresponds to a 



mixed phase. This "spheroid-halo" structure generalizes 
the "core-halo" structure of non-rotating systems. For 
even larger values of energy, the spheroid would extend 
outside the box so that the system becomes almost one- 
dimensional and depends only on the axial coordinate z 
(Fig. Eh)- Note, however, that the density contrast is ex- 
tremely low (K = p(0)/p(R) ~ 1.00022) so that the struc- 
ture is essentially homogeneous. 

3.4. Keplerian limit 

For E = E m i n , the temperature drops to zero and the sys- 
tem is equivalent to a pure polytrope with index n — 3/2. 
The structure of rapidly rotating polytropes is well-known 
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(James 1964). For rigid polytropes with index n < 0.808, 
the system is sensitive to non-axisymmetric instabilities 
of various forms. Therefore, at sufficiently high rotations, 
the spheroidal sequence bifurcates continuously to more 
complicated structures (Chandrasekhar 1969). However, 
for n > 0.808, non-axisymmetric instabilities are inhib- 
ited and the series of equilibria has a different evolution. 
The polytrope develops a cusp at the equator at a point 
where the gravitational force balances the centrifugal force 
(Fig. |2t). This is the so-called Keplerian (or break-up) 
limit corresponding to a critical angular velocity fl max 
such that Qmax-R-e = GM , where R e is the equatorial ra- 
dius. There is no equilibrium state for fl > Cl max as mass is 
ejected from the equatorial surface by the centrifugal force 
(mass shedding). The tendency to expulse some mass far 
away may be related to the formation of a Keplerian disk 
above £l ma x- In that case, we must allow for differential 
rotation. Mathematically, above £l m ax, the maximization 
problem leading to Eq. 1)11(1 has no solution that are every- 
where differentiable. The formation of an equatorial cusp 
at f2 = Umax marks the break-up of differentiability and 
the complete change of structure. For > VL max: the sys- 
tem may resemble a "cuspy condensate" + some matter 
accumulating at the box (ejected by the central body). We 
can also wonder whether non-axisymmetric structures can 
form above Q ma x- Indeed, there exists "double-cluster" 
solutions for polytropes with index n > 0.808 (Hachisu 
& Eriguchi 1984) that bifurcate dis ■continuously from the 
spheroidal sequence. On general grounds, it would be of 
interest to compare the entropy of the single cluster to 
the entropy of the double cluster. This problem is left for 
future investigations. 

For E = E m i n and O = fi max , all the mass is in the 
cuspy polytrope. As E increases along the line of maxi- 
mum rotation, the polytrope extends in size and is em- 
bedded in a disk- like halo (Fig. EF)- The maximum angu- 
lar velocity Vt max {E) decreases with E. This is because 
the size of the polytrope increases so that it is more eas- 
ily affected by rotation. At high energies, the polytrope is 
very extended and the cusp would form outside the confin- 
ing box. Therefore, the structure of the system is highly 
affected by the wall and this makes possible to achieve 
arbitrary large values of angular velocity. As a result, the 
maximum angular velocity increases with the energy when 
E is large. For high rotations and high energies, the mat- 
ter accumulates on the wall due to the centrifugal force 
(Fig. 0). The density profile is close to the law 

p = Ae^ s \ (40) 

which corresponds to the density of a rotating gas without 
self-gravity. Note that the iso-density lines are axial so that 
a complete change of symmetry has occurred from Fig.|3i 
to Fig. 0. 

3.5. Equilibrium phase diagram 

The caloric curve T(E) is plotted in Fig. for different 
values of angular velocity O. It has the typical "A-shape" 



structure previously found for non-rotating and slowly ro- 
tating systems with high cut-off. The onset of instabil- 
ity can be determined by a turning point criterion (Katz 
1978). For \l = 100, the equilibrium states are always sta- 
ble in the microcanonical ensemble (they are maxima of 
entropy at fixed mass, energy and angular momentum) but 
they are unstable in the canonical ensemble in the region 
of negative specific heats (they are saddle points of free 
energy at fixed mass, temperature and angular velocity). 
We see that, for moderate rotations, the structure of the 
caloric curve does not change dramatically with respect 
to the non-rotating case. It is just a natural continuation 
of previously reported results. However, at high rotations, 
the system is sensitive to mass shedding and there exists a 
maximum angular velocity O ma;c for each value of energy 
and temperature (Keplerian limit). 
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Fig. 3. The caloric curve of a rotating system with fi = 
100. It has a A-shape structure characteristic of self- 
gravitating systems with large cut-off (/i < hmtp = 
2670). The "gaps" in energy and temperature at high ro- 
tations are due to the absence of uniformly rotating so- 
lution above the Keplerian limit. The thin full line shows 
the £1 = case while the thick line marks the entry in the 
Kepler zone (shaded region in Figs. 0] and EJ. In between 
lines show some intermediate cases. 

The phase diagram in the (E, il) plane is plotted in 
Fig. Since the energy is assumed to be fixed, this di- 
agram corresponds to the microcanonical description . 
In continuity with the non-rotating case, the horizontal 

1 We have used the angular velocity as a control parameter 
instead of the angular momentum for commodity but also be- 
cause it is the most natural variable to characterize a rotating 
system. Indeed, the conservation of angular momentum should 
not be taken in a strict sense because it can always be satis- 
fied by ejecting a small amount of mass far away with weak 
influence on the other constraints. For similar remarks con- 
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Fig. 4. Microcanonical phase diagram in the (E, J7) plane Fig. 5. Canonical phase diagram in the (T, J7) plane for 

for [i = 100. The maximum angular velocity £l m ax(E) /i = 100. The i?-zone corresponds to an hysteretic zone 

depends on the energy. It is maximum for the T = con- where the actual phase depends on the history of the sys- 

figuration (which corresponds to an = 3/2 polytrope). tern. 
The parameters of this terminal point are Umax = 3.06 



and E„ 



1.37. 



structure of the diagram consists of three regions: a pure 
"condensed" phase for E < E conc i, a pure "gaseous" phase 
for E > E gas and a mixed "core-halo" phase for interme- 
diate energies. Vertically, the diagram proceeds with in- 
creasing flattening until the Keplerian limit is reached for 
O = ^l m ax(E). Above this line, the system is expected to 
expulse matter far away and, possibly, create a Keplerian 
disk or break into fragments. 

The phase diagram in the (T, Q) plane is plotted in 
Fig. [3] Since the temperature is assumed to be fixed, this 
diagram corresponds to the canonical description. The 
mixed phase with negative specific heats allowed in the 
microcanonical ensemble is replaced by a phase transition 
connecting the "gaseous" phase to the "condensed" phase. 
Between T c and T* the system can be in a metastable 
state. If it is initially prepared in a gaseous state, it will re- 
main gaseous until the minimum temperature T c at which 
it will collapse and become condensed. Inversely, if the 
system is initially prepared in a condensed state, it will 
remain condensed until the maximum temperature T» at 
which it will explode and become gaseous (see Fig. 0). 
Therefore, the configuration of the system in the TJ-zone 
depends on its history. This is like a hysteretic cycle in 
magnetism. 

Let us comment again on the distinction between L 
and £1. If we restrict ourselves to single-cluster axisym- 
metric solutions, as we do, we do not expect qualitative 
difference in the phase diagram according as we fix angular 

cerning the role of angular momentum in the context of vortex 
dynamics, see Brands et al. (1999). 



momentum or angular velocity. Indeed, there is a one-to- 
one correspondance between £1 and L. This is very dif- 
ferent if we allow for non-axisymmetric, or disconnected, 
solutions (such as double clusters, rings,...) that are be- 
yond the scope of the present paper. In that case, there is 
no one-to-one correspondance between £1 and L and the 
choice of the control parameter is crucial. The thermal en- 
ergy K t h = E — W — \IQ 2 is for a single central cluster 
at the same E, f2 much greater than for a double clus- 
ter with larger / and W. Thus, in (E, il) ensemble, the 
central system has more entropy (thermal random mo- 
tion) than the double cluster. Hence, the single cluster 
will overshadow the probability of the binary state (we 
thank the referee for this comment). On the other hand, 
a double cluster rotates more slowly than a single cen- 
tral cluster with same L. Therefore, in (E, L) ensemble, a 
single cluster could avoid forming an equatorial cusp by 
splitting in two pieces. Depending on L and E, the double 
cluster solution may have more entropy than the single 
one. It would be of interest to consider this problem in 
detail in the case of self-gravitating fermions. We do not 
know if "binary fermion balls" have been introduced in 
dark matter models. We can wonder whether the transi- 
tion from single to double clusters (if any) occurs at or 
before the Keplerian limit. Again, we emphasize that in 
the case of self-gravitating fermions (similar to polytropes 
n = 3/2 > 0.808), this transition would be discontinuous 
with respect to the spheroidal sequence. 
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Fig. 6. The caloric curve of a rotating system with /j, = 
10 5 . It has a Z-shape structure (dinosaur's neck) char- 
acteristic of self-gravitating systems with small cut-off 
(/•* — fJ'MTP — 2670). The ground state (T = 0, E = -138) 
is outside the range of represented energies. 



4. Caloric curves for high values of the degeneracy 
parameter (low cut-off/large system) 

We now consider the case of a high value of the degeneracy 
parameter fj, — 10 5 > hmtp = 2670 (large system) so that 
the influence of the small-scale cut-off is weaker in aver- 
age and the conditions required to set up a gravothermal 
catastrophe are fulfilled. 

The caloric curve T(E) is plotted in Fig. for differ- 
ent values of angular velocity f2. It has the typical U Z- 
shape" structure (resembling a "dinosaur's neck") pre- 
viously found 2 for non-rotating or slowly rotating sys- 



2 Of course, in a strict sense, the microcanonical caloric 
curve must be single valued everywhere because S m icro(E) = 
In g(E) and l/T micro (E) = dS m icro/dE (Gross 2001,2003). 
However, the caloric curves that we represent are richer because 
they contain saddle points SP (unstable), local maxima LEM 
(metastable) and global maxima GEM (stable) of entropy. 
Had we represented only global entropy maxima, the curves 
S(E) and f3(E) would be single-valued. They would describe 
the true statistical equilibrium states reached for t —* +oo. 
However, metastable states are of considerable importance in 
astrophysics because they correspond to the observed struc- 
tures (e.g., globular clusters) for the timescales contemplated 
(as is well known, stellar systems are not in true statistical 
equilibrium states) . Therefore, the physical caloric curves must 
take these metastable states into account and this gives rise to 
what we have nicknamed "dinosaur's necks". We recall that 
these metastable states appear only for sufficiently small cut- 
offs. For large cut-offs (see Fig. 0, there is only one global 
entropy maximum for each energy and the /3(E) curve is unival- 
ued. This is the situation considered by Votyakov et al. (2002) 
[compare their Fig. 7 to our Fig. 3]. 
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Fig. 7. Caloric curve of a self-gravitating system with 
fi = 10 5 (here Q = 0). Between E t and E c (resp. E*) the 
gaseous states (resp. condensed states) are metastable. A 
collapse occurs at E c (gravothermal catastrophe) and an 
explosion occurs at E* . By varying the energy of the sys- 
tem between E* and E c (spinodal points), we can generate 
an hysteretic cycle in the microcanonical ensemble. 

terns (Chavanis & Sommeria 1998, Chavanis 2002b, c). 
The upper branch corresponds to non-degenerate states 
(gaseous phase) and the lower branch corresponds to core- 
halo states (condensed phase). A first order microcanoni- 
cal phase transition is expected to occur at E t (//) since at 
this energy the two phases have the same entropy (Fig. EJ ■ 
This gravitational phase transition is marked by a discon- 



gas 



tinuity of temperature and specific heats. If E t > E, 
(where E gas is the first turning point of temperature), the 
specific heats passes from positive to negative values. If 
Et < E gas , the specific heats is always negative at the 
transition (the crossing point occurs for /i ~ 1.10 10 4 ). 
However, this first order phase transition may not take 
place in practice because the gaseous states with energy 
E c < E < E t are metastable (local entropy maxima) and 
long-lived (Katz & Okamoto 2000, Chavanis & Ispolatov 
2002). The collapse (gravothermal catastrophe) will set in 
at, or near, the critical energy E c (Antonov energy) at 
which the gaseous phase disapears. This corresponds to 
a zeroth order microcanonical phase transition accompa- 
nied by a discontinuity of temperature and entropy. The 
collapse stops when the central region becomes degener- 
ate and forms a dense nucleus. The resulting configuration 
has a "core-halo" structure. Contrary to the collapse in the 
canonical ensemble, the nucleus contains only a moderate 
fraction of mass. In the classical limit fi — > +oo, it reduces 
to a single binary (see Sec.EJ). The rest of the mass is di- 
luted in a hot envelop (with almost uniform density) that 
is held by the walls of the box. For an open system, the 
halo would be dispersed at infinity so that only the de- 
generate nucleus would remain. This behaviour, with the 
expulsion of an envelope, is more consistent with the pro- 
cess of white dwarf formation than a complete collapse 
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Fig. 8. Microcanonical phase diagram for p = 10 5 . The H- 
zone corresponds to an hysteretic zone where the actual 
phase depends on the history of the system. 

of the system at fixed T. This is an astrophysical evi- 
dence that the microcanonical ensemble is more appro- 
priate than the canonical one. The condensed states with 
energy E t < E < E* are also metastable and long-lived. 
An "explosion" , reversed to the gravothermal catastrophe, 
will occur above E* if the system is initially prepared in 
the condensed phase (see Fig. EJ- By varying the energy 
between E* and E c (microcanonical spinodal points) we 
can generate an hysteretic cycle in the microcanonical en- 
semble (Chavanis 2003b). 

The phase diagram describing the microcanonical 
phase transition is represented in Fig. |S1 It can be com- 
pared with the phase diagram describing the canonical 
phase transition for p = 100 (see Fig. [2J , the role of energy 
E and temperature T being reversed. In particular, there 
exists a ii-zone in which the structure of the system (gas 
or condenstate) depends on its history For smaller ener- 
gies (E < E c ), the system is in a mixed phase with nega- 
tive specific heats. For even lower energies E con d ~ —117 
(not represented on the diagram), we leave the region of 
negative specific heats and enter the pure condensed phase 
as in Fig. 0| More details on these phase transitions can be 
found in Chavanis (2002b). It is shown in particular that 
the microcanonical first order phase transition disapears 
for p < pmtp — 2670. 

5. Gravothemal catastrophe for classical 
self-gravitating particles (no cut-off) 

We now consider the limit p = +oo (or h = 0) correspond- 
ing to classical particles. In that limit, the particles do not 
feel any kind of exclusion or small-scale cut-off. This is the 
situation relevant for globular clusters for which the size 
of the stars does not matter. 



In the non degenerate limit, the caloric curve has a 
classical spiral shape (Fig. |5J). We recall that, in the ab- 
sence of small-scale cut-off, the equilibrium states are 
only metastable as they correspond to local maxima of 
the thermodynamical potential (entropy or free energy). 
These metastable states play however an important role 
as they correspond to the observed structures (see dis- 
cussion in Chavanis 2003a) . In the microcanonical ensem- 
ble, the gravothermal catastrophe sets in at the Antonov 
energy E c . The microcanonical phase diagram in (E, f2) 
plane is drawn in Fig. and shows in particular the evo- 
lution of the Antonov critical energy E c (£l) with the rota- 
tion. The critical energy E C (SY) increases with rotation so 
that gravitational collapse occurs sooner than in the non- 
rotating case. This confirms the results of the perturbative 
approach (Chavanis 2002c). Qualitatively, the critical en- 
ergy increases due to the rotational term \IVL 2 > in 
the expression of E (see Sec. 12.50 . The result of the col- 
lapse is to form a binary star with large binding energy. 
Since total energy is conserved, the gravitational energy 
released by the binary is redistributed in a dilute halo 
in the form of thermal energy. This evolution is natural 
in the microcanonical ensemble since it can produce very 
high (ideally diverging) values of entropy due to the di- 
vergence of the temperature. It is also vindicated by nu- 
merical simulations of the orbit averaged Fokker-Planck 
equation (Cohn 1980) or fluid equations (Lynden-Bell & 
Eggleton 1980). These simulations show that the collapse 
of globular clusters is self-similar and leads to a finite time 
singularity with a density profile p ~ r~ a with a ~ 2.2. 
The mass contained in the core goes to zero (indicating 
that the collapse is halted by the formation of binaries 
with mass 2m <C M) while the central temperature raises 
to +oo since a > 2. Binary- like structures are also found 
in the microcanonical model of self-gravitating Brownian 
particles (Sire & Chavanis 2003). 

In the canonical ensemble, the collapse sets in at 
the critical temperature T c . The evolution of this criti- 
cal temperature with the angular velocity is represented 
in Fig.r^| It is found that T c (f2) decreases with f2 so that 
gravitational collapse is delayed in the canonical ensemble. 
This corroborates the results found with the perturbative 
approach (Chavanis 2002c). The result of the collapse is 
to form a Dirac peak, i.e. a state where all the mass is con- 
centrated at the center ("black hole"). This evolution is 
natural in the canonical ensemble since it can produce very 
high (in fact diverging) values of free energy J = S — j3E 
due to the divergence of the potential energy. This isother- 
mal collapse can be followed dynamically by solving the 
Euler-Jeans equations for a gaseous system with an equa- 
tion of state p = pT where T is fixed (Penston 1969). It 
can also be studied analytically in the canonical model of 
self-gravitating Brownian particles (Chavanis et al. 2002). 
In this model, there is first a finite time singularity lead- 
ing to a r~ 2 density profile. Then, the collapse continues 
after the singularity is formed and the Dirac peak appears 
in the post-collapse regime due to the accretion of mass 
(Sire & Chavanis 2002,2003). We expect this two-stage 
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Fig. 10. Dependence of the Antonov critical energy and critical temperature with rotation. For high energies or 
high temperatures, the system is in a gaseous (metastable) phase. Below a critical energy or temperature, it enters 
in the collapsed phase. In the microcanonical ensemble (fixed E), the collapse leads to a single binary (or a small 
group of stars) surrounded by a hot dilute halo. In the canonical ensemble (fixed T), the collapse leads to a Dirac 
peak containing all the mass ("black hole"). This Dirac peak forms in the post- collapse regime of the dynamical 
evolution. For small rotations, the critical energy and critical temperature behave as E c (fl) = E c (0) + 0.175f2 2 and 
1/T C (Q) = 1/T C (0) + 0.552H 2 (dotted lines). This scaling is valid typically for Q < 0.5. 



behaviour to persist in the case of more relevant dynami- 
cal models. This may be an important point for numerical 
simulations in astrophysics and cosmology because if we 
only consider the first regime we conclude that the core 
contains no mass (while its density is high and scales as 
r~ 2 ). However, in a post-collapse regime, a Dirac peak 
is formed which contains a macroscopic mass at a single 
point. We can wonder whether this mechanism (appar- 
ently not reported before) is related to the existence of 
black holes at the center of galaxies. In reality, small-scale 
cut-offs must be introduced (e.g., quantum degeneracy) 
and the black hole singularity is replaced by a "fermion 
ball" with a high mass. Note that in the microcanonical 
ensemble, there is no black hole singularity for classical 
point masses (the Dirac peak is replaced by a binary). 
Interestingly, for self-gravitating fermions, the condensate 
(fermion ball) contains a moderate but macroscopic frac- 
tion of mass (Chavanis & Sommeria 1998). This may ex- 
plain the formation of massive objects at the center of 
galaxies in a microcanonical framework. The microcanon- 
ical description is the most relevant as it accounts for the 
transfers of energy between the core and the halo contrary 
to isothermal models. 

The divergence of entropy in the microcanonical en- 
semble due to the formation of a binary embedded in a hot 
halo is implicit in the work of Antonov (1962). It is also 
clear from the arguments given by Padmanabhan (1990) 



to show that the density of state diverges as we approach 
two particles to each other provided that thermal energy 
can be redistributed in a halo (this demands N > 3). 
The divergence of free energy in the canonical ensemble 
by forming a Dirac peak has been shown rigorously by 
Kiessling (1989) and heuristically by Chavanis (2002e). 
The inequivalence of statistical ensembles regarding the 
formation of binaries or Dirac peaks is further discussed 
in Appendices A and B of Sire & Chavanis (2002). In 
this approach, a family of distributions functions is con- 
structed which provokes the divergence of entropy or free 
energy when a physical parameter is varied. These results 
can also be obtained by considering the small cut-off limit 
fj, — > +oo in the analytical model of phase transitions for 
self-gravitating systems proposed by Chavanis (2002b). 

6. Phase diagram of non-rotating systems 

In previous sections, we have presented the equilibrium 
phase diagram of self-gravitating fermions as a function 
of angular velocity f2 for a fixed value of degeneracy pa- 
rameter (j,. We shall now discuss the equilibrium phase 
diagram of self-gravitating fermions (or hard sphere sys- 
tems) as a function of /j, for Q = 0. These phase diagrams 
have not been published before. 

The deformation of the caloric curve T(E) as a func- 
tion of the degeneracy parameter \i is described in detail 
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Fig. 9. Caloric curves corresponding to classical point 
masses (fi = +oo) for different values of angular velocity. 
Below E c , the gravothermal catastrophe leads to the for- 
mation of a binary star (microcanonical ensemble) . Below 
T c , the isothermal collapse generates a Dirac peak (canon- 
ical ensemble). 



by Chavanis (2002b). This work completes earlier studies 
on phase transitions in self-gravitating systems by showing 
the small cut-off limit, characterized by an unwound spiral 
("dinosaur's neck"), leading to a first order microcanoni- 
cal phase transition (earlier works essentially described the 
canonical first order phase transition). Typical curves il- 
lustrating first order canonical and first order microcanon- 
ical phase transitions are shown in Figs.n an d|3 The equi- 
librium phase diagram of self-gravitating fermions can be 
directly deduced from these results by identifying charac- 
teristic energies and temperatures. In the canonical ensem- 
ble, we note T t the temperature of transition (determined 
by the equality of the free energies of the two phases), 
T c the end point of the metastable gaseous phase (first 
turning point of temperature) and T* the end point of the 
metastable condensed phase (last turning point of temper- 
ature). The temperature of transition T t is determined by 
a horizontal Maxwell construction (Chavanis 2003b). The 
canonical phase diagram is shown in Fig. 1111 It shows in 
particular the canonical critical point \iqtp = 83 at which 
the canonical first order phase transition disappears and 
is replaced by a localized second order phase transition. 

In the microcanonical ensemble, we note E t the en- 
ergy of transition (determined by the equality of the en- 
tropy of the two phases), E c (Antonov energy) the end 
point of the metastable gaseous phase (first turning point 
of energy) and E* the end point of the metastable con- 
densed phase (last turning point of energy). The energy 
of transition E t is determined by a less standard verti- 
cal Maxwell construction (Chavanis 2003b). We also de- 
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Fig. 11. Canonical phase diagram in (T, fi) plane for £1 = 
0. The H-zone corresponds to an hysteretic zone where 
the actual phase depends on the history of the system. 

note by E gas the energy at which we enter in the mixed 
phase with negative specific heat (first turning point of 
temperature) and E conc [ the energy at which we leave the 
mixed phase and enter the condensed phase (last turn- 
ing point of temperature). We also introduce the min- 
imum energy E m i n (ground state). The microcanonical 
phase diagram is shown in Fig. 1121 It shows in particular 
the microcanonical critical point \imtp = 2670 at which 
the microcanonical first order phase transition disappears. 
The structure of the equilibrium phase diagrams can be 
easily understood in the light of the preceding discussion 
(see Chavanis 2002b for more details). We note that the 
microcanonical phase diagram is more complex than the 
canonical one due to the existence of a negative specific 
heat region. We recall that canonical stability (maximum 
of J at fixed M, /3) is a sufficient, albeit not necessary, con- 
dition of microcanonical stability (maximum of S at fixed 
M, E). Since canonical equilibria are always realized as 
microcanonical equilibria, they constitute a sub-domain 
of the microcanonical phase diagram. Hence, the micro- 
canonical ensemble is richer than the canonical one, as is 
well-known. 

7. Conclusion 

We have calculated single-cluster axisymmetric statisti- 
cal equilibrium states of self-gravitating fermions for ar- 
bitrary values of rotation and temperature (or energy). 
The problems associated with the divergence of the grav- 
itational potential as r — > have been avoided by con- 
sidering quantum effects. This small-scale regularization 
is a physical one because the Pauli exclusion principle is 
a fundamental concept and fermionic matter is relevant 
in astrophysics. Isothermal and polytropic configurations 
(with index n = 3/2) are obtained as particular limits of 
our model (T — ► +oo and T = 0). On the other hand, de- 
pending on the importance of rotation, we have obtained 
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Fig. 12. Microcanonical phase diagram in (E, p) plane for 
O = 0. The phase diagram in MCE is more complex than 
in CE due to the existence of the negative specific heat 
region that is forbidden in CE. 



different types of structures such as fermion balls, fermion 
spheroids and cuspy structures. In the mixed phase, these 
condensed objects are surrounded by a diffuse isothermal 
halo. Clearly, the main drawback of our study is the ne- 
cessity to enclose the system within an artificial box so 
as to prevent its complete evaporation (which is the true 
statistical equilibrium state for self-gravitating systems). 
A future extension of our study is to consider more real- 
istic truncated models in order to avoid the artifice of a 
material box and take into account incomplete relaxation 
or tidal effects. We note, in passing, that the distribution 
function / = /(ej) predicted by statistical mechanics can- 
not account for the triaxial structure of elliptical galax- 
ies since no bifurcation from axisymmetric (spheroidal) 
to non-axisymmetric (ellipsoidal) structures are expected 
for uniformly rotating compressible bodies. The distribu- 
tion function must depend on isolated integrals different 
from the single Jacobi integral (Contopoulos 1960). This 
is clearly a manifestation of incomplete relaxation but the 
construction of such isolated integrals is a highly compli- 
cated task and demands a very delicate analysis of the 
process of formation of elliptical galaxies. 

We emphasize that the nature of the small-scale cut- 
off is of considerable importance for rotating objects, re- 
garding the structure of the condensed phase, because 
it can generate or inhibit non-axisymmetric instabilities. 
For example, the ordinary Fermi-Dirac distribution used 
in this paper and the model of Votyakov et al. (2002) 
have the same high temperature limit (an isothermal gas 
with p = pT) but a different low temperature limit. The 
Fermi gas at T = corresponds to the equation of state 
p = Kp 5 ^ 3 (polytropic core). The equation of state cor- 
responding to a Fermi-Dirac distribution in configuration 
space p = p mQ;E /[l + Acxp(/3$ e// )] isp = Tln(l- pj ' p max ), 
as can be obtained from the condition of hydrostatic equi- 



librium 3 . For T = 0, p is a step function with p — > p max 
for r < i?» (solid core) and p = otherwise. In the 
former case (fermions), there is no bifurcation to non- 
axisymmetric structures at high rotations in continuity 
with the spheroidal sequence. Instead, the system devel- 
ops an equatorial cusp (the possibility of discontinuous 
bifurcations to non-axisymmetric structures will be con- 
sidered in a future work). This cusp does not exist for 
homogeneous bodies (or hard sphere systems) and non- 
axisymmetric bifurcations occur continuously at high ro- 
tations. It is surprising that the equation of state corre- 
sponding to the model of Votyakov et al. (2002) does not 
coincide with the Van dcr Waals equation of state, as in 
the work of Aronson & Hansen (1972), since these authors 
treat their system as a hard sphere gas 4 . It is also surpris- 
ing that they do not obtain a richer variety of structures in 
the low temperature/energy regime (pear shaped, dumb- 
bell,...) since rotating homogeneous bodies possess a rich 
bifurcation diagram as investigated since the 19-th cen- 
tury (Chandrasekhar 1969). These remarks show that the 
small-scale regularization must be given precise consider- 
ation depending on the situation contemplated. 

The self-gravitating Fermi gas has several astrophysi- 
cal applications: (i) Lynden-Bell's type of degeneracy may 
be relevant for galactic nuclei and dark matter (Lynden- 
Bell 1967, Kull et al. 1996, Chavanis & Sommcria 1998, 
Stiavelli 1998, Leeuwin & Athanassoula 2000). (ii) Pauli's 
exclusion principle is relevant for white dwarfs and neu- 
tron stars (Chandrasekhar 1942, Hertcl & Thirring 1971). 
In that context, the Fermi gas at non-zero temperature 
provides a simple theoretical model describing a phase 
transition from a "gaseous star" (for T > T c ) to a 
white dwarf or a neutron star (for T < T c ). To obtain 
a better description of these objects, special relativity 
must be taken into account for massive white dwarf stars 

3 A Fermi-Dirac distribution in configuration space also 
arises in the statistical mechanics of two-dimensional vortices 
(see Chavanis 2002a for a review on the analogy between stellar 
systems and 2D vortices). In that respect, it may be relevant to 
mention the work of Chen & Cross (1996) who found "double- 
vortex" equilibrium solutions in a circular domain when the 
conservation of angular momentum is accounted for. This can 
be viewed as the 2D hydrodynamical version of the "double- 
star" structure found by Votyakov et al. (2002). Note that the 
presence of a confining box is crucial to maintain the double- 
vortex solution. Similarly, for self-gravitating systems at non- 
zero temperature, an artificial box is needed to maintain the 
double-star structure. 

4 In fact, Votyakov et al (2002) prevent overlapping of parti- 
cles in physical space, as in a lattice gas model. This is a sim- 
ple way to limit the local density p(r) of the self-gravitating 
gas. Hence, their model is expected to represent stars at densi- 
ties up to the ignition of their hydrogen burning, when further 
gravitational collapse is halted for a while (Gross 2003). In 
this context, it may be relevant to consider a large small-scale 
cut-off, as they do. Note that their model does not correspond 
to the Lynden-Bell (1967) statistics, which is a Fermi-Dirac 
distribution in phase space appropriate to collisionless stellar 
systems (e.g., Chavanis 2002d). 
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(Chandrasekhar 1931) and general relativity is required 
in the case of neutron stars (Oppenheimer & Volkoff 
1939). (iii) It has been proposed recently that dark mat- 
ter could be made of a collisionless gas of massive neutri- 
nos. Therefore, the relevant distribution function is again 
the Fermi-Dirac distribution (see, e.g., Chavanis 2002f). 
The degeneracy is due either to the Liouville theorem (in 
the context of violent relaxation) or to Pauli's exclusion 
principle (if quantum effects are relevant) . By cooling be- 
low a critical temperature, this neutrino gas is expected 
to undergo a phase transition to a compact object called 
a "fermion ball" . It has been proposed that fermion balls 
could provide an alternative to black holes at the center of 
galaxies (Bilic & Viollier 1997). (iv) Incidentally, our study 
can also be connected to models of rotating stars (e.g., 
Roxburgh et al. 1965). Low mass stars (M < 1.5 M ) of 
the main sequence have a radiative core and a convective 
envelope. Massive stars (M > 1.8 M Q ) have the opposite 
configuration. Now, the convective region corresponds to 
a polytrope of index n = 3/2. Coincidentally, this case is 
precisely covered by our study. 

On the other hand, hard spheres models have less di- 
rect applications in astrophysics because the inter-particle 
distance is always considerably larger than their size. 
Indeed, globular clusters can be either in dilute metastable 
equilibrium states (Michie 1963, King 1966) or in a col- 
lapsing phase ending up with the formation of binaries 
(Henon 1961). These binaries can release sufficient energy 
to drive a re-expansion of the system (Inagaki & Lynden- 
Bell 1983, Bettwieser &; Sugimoto 1984) so that complete 
collapse is prevented in practice. In each case, the size 
of the stars does not matter. The size of the atoms in 
stars does not matter neither because radiative transfers 
keep the star in a gaseous phase far from condensation. 
When a star collapses (because it has burned its nuclear 
fuel) , gravity is balanced by quantum pressure not by the 
pressure of a hard sphere gas. In that respect, it may be 
useful to recall that the hard sphere model of Aronson & 
Hansen (1972) was presented as a crude model for neutron 
stars where the hard spheres were introduced to mimick 
quantum degeneracy. The hard-sphere model can however 
find physical applications in the context of planet forma- 
tion resulting from the gravitational collapse of dust parti- 
cles in the solar nebula (Chavanis 2000). It has also been 
used by Stahl et al. (1994) to describe the formation of 
stars or planetoids (considered as the solid phase) result- 
ing from the condensation of a molecular cloud (consid- 
ered as the gaseous phase). In this context, the explosion, 
reverse to the collapse, may be related to the supernova 
phenomenon. 

Finally, our study has a direct relevance in terms of 
fundamental statistical mechanics. Indeed, the different 
types of phase transitions arising in the self-gravitating 
Fermi gas (Chavanis 2002b) display a very rich structure: 
inequivalence of statistical ensembles, negative specific 
heats, convex intruder in the entropy-energy curve, canon- 
ical and microcanonical first order phase transitions, sec- 
ond order phase transition, zeroth order phase transition, 



canonical and microcanonical critical points, metastable 
states, hysteretic cycles,... Such properties are probably 
common to other systems with long-range interactions. In 
fact, the self-gravitating Fermi gas forms a generic exam- 
ple of systems displaying codimensions and 1 singular- 
ities in the classification of phase transitions (for long- 
range interactions) recently proposed by Bouchet & Barre 
(2003). Therefore, the detailed study of phase transitions 
in the self-gravitating Fermi gas at non-zero temperature 
is of interest both for statistical mechanics and astro- 
physics. This point is important because the statistical 
mechanics of self-gravitating systems was considered until 
recently, as a curiosity, not to say a fallacy, by statistical 
mechanicians. 

An extension of our study would be to compute more 
complex configurations of rapidly rotating self-gravitating 
fermions and n — 3/2 polytropes (double clusters, rings 
and core-ring sequences) in both microcanonical (fixed E, 
M and L) and canonical (fixed T, M and S~2) ensembles. 
These solutions are expected to bifurcate discontinuously 
from the spheroidal sequence in contrast to hard sphere 
systems and homogeneous bodies (n — 0). It would be im- 
portant (but difficult) to investigate the stability of these 
various configurations. For rigid bodies {n < 0.808 poly- 
tropes and hard sphere systems at low temperature), the 
spheroidal sequence becomes unstable (saddle points of 
S or J) above a certain rotation at which the branch 
of non-axisymmetric structures appears. For high rota- 
tions, double clusters are global maxima of S or J and 
the ring sequence is always unstable (saddle points of S 
or J). For compressible bodies (n > 0.808 polytropes and 
fermionic matter at low temperature), the spheroidal se- 
quence is expected to remain stable (maxima of S or J) 
until the Keplerian limit since no continuous bifurcations 
exist. However, it is possible that, for high rotations, the 
single cluster structures computed in this paper are only 
metastable (local maxima of S or J) and that the true 
equilibrium states (global maxima of S or J) are double 
clusters or core- ring configurations. In case when several 
stable solutions exist for the same values of E, L (or j3, 
ft), it may be important to compare their entropy and 
free energy. However, the selection of the physical equilib- 
rium state will depend more on a notion of basin of attrac- 
tion (see, e.g., Chavanis et al. 2002 in a related context) 
than whether it is a global or a local maximum of S or 
J. Thus, a dynamical study is required. We hope to come 
to these problems in future communications. The present 
paper is only a first step in that direction. To be more re- 
alistic, the box could be avoided by using the distribution 
function l|15(l which provides a more physical confinement 
of the system. We also emphasize that in a large domain, 
the conservation of angular momentum can be satisfied by 
ejecting a weak amount of mass at large distances, for ex- 
ample in the form of a Keplerian disk. This is the typical 
structure of a protoplanetary nebula. Indeed, during the 
formation of a (single) star by gravitational collapse of a 
molecular cloud, about 99% of the intial angular momen- 
tum is spread in a protoplanetary disk while 99% of the 
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mass remains in the star itself whose internal structure is 
hardly affected by rotation. 
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